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It has been shown recently that Dirac operators satisfying the Ginsparg- Wilson relation provide a solution of 
the chirality problem in QCD at finite lattice spacing. We discuss different ways to construct these operators and 
their properties. The possibility to define lattice chiral gauge theories is briefly discussed as well. 



1. Introduction 

Chiral symmetry plays an essential role in par- 
ticle physics. In QCD it is a global symmetry 
which forbids an additive quark mass renormal- 
ization and whose spontaneous breakdown pro- 
vides us with (nearly) Goldstone bosons with 
their specific interactions. 

On the other hand, in electroweak interactions 
chiral symmetry plays an even more fundamental 
role: it is a local gauged symmetry. Its presence 
is necessary to provide a renormalizable theory. 

Accordingly, when trying to regularize these 
theories on the lattice (which is the only non- 
perturbative regularization known for gauge the- 
ories) one meets two basic problems. As it is (or 
rather was) well known one cannot realize chiral 
symmetry exactly on the lattice without sacrific- 
ing some even more important principles, like lo- 
cality or absence of extra physical particles. This 
is the content of the famous no-go theorem by 
Nielsen and Ninomiya 0. In Wilson's formula- 
tion of lattice QCD chiral symmetry is explicitly 
broken by an irrelevant term which does not affect 
the continuum limit. It causes, however, 0(a) lat- 
tice artifacts and other unwanted effects at finite 
lattice spacing. 

For the electroweak theory, on the other side, 
the problem of chiral symmetry on the lattice is 
a principal one: is it possible to define a chiral 
gauge theory non-perturbatively? This second 
problem is much more difficult, and assumes that 
the first one - a solution for global chirality - has 
been successfully solved. 
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Several ways to solve these problems have been 
proposed. These were mostly covered on pre- 
vious lattice conferences. Here I mention only 
the domain wall fermions ||J| and the related 
overlap formalism 0J|. On can state that these 
approaches provide a solution to the problem of 
global chirality. The numerical results using do- 
main wall fermions Q are encouraging. In the 
overlap formalism also important steps were made 
towards the solution of the gauged chiral symme- 
try (for references see Ref. ||). 

It was unclear, however, to what extent are 
the theoretical constructions used in these ap- 
proaches necessary, or what is the basic ingredient 
behind the chiral properties observed. As it be- 
came clear recently, the basic relation which has 
to be satisfied to realize chiral symmetry at finite 
lattice spacing, is the Ginsparg- Wilson (GW) re- 
lation H suggested a long time ago! 

In this talk I shall concentrate on the role of 
the GW relation and its consequences. Some of 
the results were obtained earlier in the approaches 
mentioned above, but even these are more trans- 
parent within the new framework. 

1.1. The Nielsen-Ninomiya no-go theorem 

Consider the free fermionic action 

S F = a i Y / ^)D(x-yMy), (1) 

describing a massless fermion on the lattice. 
(Colour, Dirac and flavour indices are sup- 
pressed.) 

The desirable properties of the Dirac operator 
D should be 

(a) D(x) is local (bounded by Ce -7 ^') 
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(b) D(p) = ij^Pft + 0{ap 2 ) for p < ir/a 

(c) D(p) is invertible for p / (no massless 
doublers) 

(d) j 5 D + D l5 = 

Here D(p) is the Fourier transform of D(x). 

The theorem by Nielsen and Ninomiya Q 
states that properties (a)-(d) cannot hold simul- 
taneously. 

Following Wilson, to avoid the massless dou- 
blers present in the naive discretization of the 
Dirac operator, one introduces a chiral symmetry 
breaking, irrelevant operator. The corresponding 
Wilson-Dirac operator is given by 

A^^v^ + vp-iav;;^, (2) 

where V M (V*) is the forward (backward) lattice 
derivative. 

Due to the explicit breaking of chiral invari- 
ance, however, some extra problems appear: 

(1) Chiral symmetry is broken at finite lattice 
spacing by 0(a) lattice artifacts and it is recov- 
ered only in the continuum limit. 

(2) One has to fine tune the bare quark mass 
since there is an additive quark mass renormal- 
ization. 

(3) Operators of different chiral representations 
get mixed, etc. 

1.2. Locality 

One of the conditions of the no-go theorem 
was locality. The natural definition (used here) 
is that the Dirac operator in an arbitrary gauge 
field background falls off exponentially i.e. satis- 
fies the bound 

\\D(x,y;U)\\<Ce-^-y\, (3) 

where C and 7 > are independent of the gauge 
field U. In addition, it is also assumed that 
D(x,y;U) depends negligibly on U's far apart 
from x and y, i.e. that 



also falls exponentially in \x — z\, \y — z\ and \x — 
y\. It is generally assumed that such interactions 
satisfy universality. 



Note that it is physically too restrictive to con- 
sider only nearest neighbour interactions or even 
those with strictly finite support. In a real sys- 
tem the interaction coefficients die away expo- 
nentially rather then having a finite number of 
strictly nonzero terms. 

On the other hand, interactions decreasing 
more slowly then exponentially should be consid- 
ered as non-local and non-acceptable, since uni- 
versality is not expected to hold for these. Some 
non-local interaction can give the same correla- 
tions as a local one (e.g. when a bad RG trans- 
formation is used, or when some physical fields 
are integrated out). This is, however, no suffi- 
cient excuse to use them - a slight change in the 
interaction could modify the picture drastically. 

Theoretically any decay rate 7 in eq. (||) is al- 
lowed. In practice, however, one works at finite 
lattice spacing and 7 > ma must hold. For that 
reason it is important to have a sufficiently small 
interaction range. 

In principle it could happen that the decay rate 
of D{x, y\ U) is not uniformly bounded by an ex- 
ponential, but there exists a bound j(U) (with 
inf[/7({7) = 0). In general, such interactions 
should be considered non-local. However, it may 
happen that those configurations when "/(U) < 70 
are so strongly suppressed that they practically 
'never' occur in a MC simulation. This is a po- 
tentially dangerous possibility, but could perhaps 
be accepted with extra caution. 

2. The Ginsparg- Wilson relation 

Back in 1982 Ginsparg and Wilson || sug- 
gested a way to avoid the no-go theorem and 
preserve consequences of the chiral symmetry. 
They suggested to require instead of relation 
75!? + -D75 = the following milder condition 

75^ 1 + Z?- 1 7 5 =a2i? 75 . (5) 

Here a is the lattice spacing and R is a local oper- 
ator. For simplicity it will be assumed also that 
R is trivial in Dirac indices. That R commutes 
with 75 follows from eq. ([|). The locality of R 
expresses the requirement that the propagating 
states are effectively chiral, so at distances larger 
than the range of R its presence is not felt. Ob- 
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viously, this is a highly nontrivial condition since 
the propagator D~ x on the lhs. is non-local. (The 
coefficient 2 on the rhs. is for historical reasons 
only.) 

Accordingly, the Dirac operator D should sat- 
isfy the Ginsparg- Wilson relation 

75-D + D To = aD2Rj 5 D (6) 

In Ref. H this relation originates from RG con- 
siderations applied to the free fermionic action 
and has been generalized to fermions in a gauge 
field. The local operator R comes from a chiral 
symmetry breaking term in the block transforma- 
tion only, not from the original action of the RG 
procedure. For that reason the physics described 
by D is expected to be chirally invariant. (For 
a chirally symmetric block transformation D be- 
comes chirally invariant but non-local, in accor- 
dance with the no-go theorem.) 

Two comments are in order here: 

(1) The naive commutation property for D is 
recovered in the continuum limit a — > 0. 

(2) The rhs. of eq. @ is zero on solutions, i.e. 
for Dip = 0. 

Ginsparg and Wilson suggested that a solution 
to eq. (^) is the mildest way to break chiral in- 
variance: it should maintain the physical conse- 
quences of chiral invariance (soft pion theorem, 
etc.). However, no solution was found for the in- 
teracting case (QCD), and the paper has been 
practically forgotten for 15 years! 

Last year Peter Hasenfratz realized [Q that the 
fixed point (FP) action (or classically perfect ac- 
tion) for QCD § satisfies the GW relation. This 
observation revived the interest in this relation. 

2.1. The fermionic FP action 

Let me recall the definition of the FP action. 
For the pure gauge part of the action one has the 
saddle point (i.e. classical) equation 

S F g p (V) = min{S g FP ([/) +T g (V,U)} . (7) 

Here the gauge field V lives on the coarse lat- 
tice, U on the fine one, and T g (V, U) is the block- 
ing kernel. This is a recursive definition, but its 
recursive evaluation converges rapidly since the 
minimizing field U — U(V) is much smoother 
than the original V. 



In the classical approximation the fermionic 
part remains quadratic, and the corresponding 
propagator satisfies the recursion relation 

D^(n,n';V) = -5 nn , (8) 

K 

+ ^2 ^(x, n; U)Dpp(x, y; U)u(y, n'\ U) 

Here n, n' label the sites on the coarse lattice, x, y 
on the fine lattice, while on the rhs. U — U(V), 
the minimizing configuration of the pure gauge 
problem is taken. The term proportional to 1/k 
comes from the chiral breaking term in the block- 
ing kernel for fermions; Lu(x,n;U) characterizes 
the averaging procedure for the fermions in RG 
transformation - it is assumed to be trivial in 
Dirac indices and usually it only extends over the 
hypercube. Applying eq. (|J) recursively and us- 
ing the fact that the original D (in the continuum) 
anti-commutes with 75, it is easy to see that -Dpp 
satisfies eq. (||) with 

R(U) = - (local operator) . (9) 

K 

(In fact, for the RG transformation used, R nn i cx 
5 nn i or it lives on the hypercube, i.e. is even ultra- 
local.) 

In Ref. H Hasenfratz has demonstrated that, 
indeed, the GW relation implies the chiral Ward 
identities. The argument is the following: When 
the symmetry breaking part appears in some ex- 
pressions (integrated out in fermionic fields) 

(...^DR lb Dip...) (10) 

then the D~ 1 factors from contractions of the 
fermionic fields are cancelled by the D's on both 
sides in the breaking term, leaving only the local 
operator R. As a consequence, the symmetry re- 
lation d^J^f^ ~ is satisfied up to a contact term 
- i.e. one obtains the Ward identities. In this 
paper it has also been shown that the proper or- 
der parameter for spontaneous chiral symmetry 
breaking is given by 

(^xV>x) sub = ($J>x ~ Rxx) ■ (11) 

It has been proven as well that the current is not 
renormalizcd (Za = Zy = 1) and operators in 
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different chiral representations do not mix. I shall 
discuss these issues later, in the light of newer 
developments. 

In Ref. jlO| it has been shown that there is an 
exact Atiyah-Singer index theorem on the lattice, 
and the fermionic spectrum (in arbitrary gauge 
background) has nice chiral properties. The dis- 
cussion was based on the FP action, but it used 
mostly only the GW relation to establish these 
facts. 

As mentioned in the introduction, the domain 
wall and the related overlap formalism also pro- 
duced Dirac operators with nice chiral properties. 
Neuberger has pointed out that the Dirac op- 
erator obtained earlier |12| also satisfies the GW 
relation. To be on more general ground, I discuss 
first the overlap Dirac operator. 

2.2. Neuberger's Dirac operator 

This is an explicit construction. Introduce fl^| 

A = l-aD w , (12) 

where Z? w is the standard massless Wilson-Dirac 
operator given by eq. (^). Note that 

75-Dw75 = Dl, . 
Then define 



(13) 



a 



It satisfies the GW relation 
75 D + D75 = aDj 5 D 



(14) 



(15) 



(i.e. eq. (||) with R = 1/2; sometimes this spe- 
cial case will be referred to as the GW relation). 
Introducing V through 

V=l-aD, (16) 

together with the relation 75-D75 = one con- 
cludes that eq. (15) is equivalent to 



V*V = l, 

i.e. V is unitary. In other words, 
1 



D 



-(1-V) 



(17) 



(18) 



with a unitary V defines a general solution to 
eq. ([l5|). Of course, eq. ( |l8| ) in itself is not enough 



to define an acceptable Dirac operator. For the 
free case the overlap Dirac operator has the prop- 
erties (a)-(c) listed in the no-go theorem. 

It will be argued later that D remains local in a 
gauge field background as well. As a consequence, 
D(x, y, U) defines an acceptable lattice Dirac op- 
erator (properties (a)-(c) of no-go theorem). 

Let me make a few simple comments. (Some 
of these straightforward observations in this or 
modified form have been made independently by 
several people and I shall not always give a refer- 
ence. See, however, Refs. |l3|.14].) 

(1) Instead of D w in eq. (12) one could start 
with any acceptable Dirac operator Dq (massless, 
no doublers, local). 

(2) If the original Dq satisfies the GW relation 
eq. (|l5|) then A^A = 1 hence D = D , i.e. the 
operator reproduces itself. 

(3) H = 75A = 75(1 — -D w ) is hermitian and 
A^A = H 2 . As a consequence 



D = ± (1 - l5 e(H)) 



(19) 



where e(x) is the sign function. 

(4) For A = A(fi) = [i-aD w (where < \i < 2) 
one has 

1 ( ,1 
D = n-[l- A—=_ 

a V \[ATAj 
which satisfies a slightly modified relation 
1 



(20) 



75 D + D75 



I' 



-aD-f 5 D 



(21) 



The choice of fi 7^ 1 will be useful as discussed 
later. Note that for /i < there are no mass- 
less fermions at all, while for /x > 2 there are too 
many, since the modes around X(D W ) — 2 (the 
massive doublers) also become massless modes of 
D. Therefore only < /1 < 2 produces an accept- 
able Dirac operator. 

(5) It is also easy to generalize to the case of 
arbitrary local R: 



A = l 
and 



'2RD^2R, 



D = - 

a 



1 1 



1-A- 



(22) 



(23) 



2R V s/ATA) \/2R 

This satisfies eq. (|^). For simplicity, however, we 
shall consider mostly the case R = 1/2. 
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3. The index theorem on the lattice 

For a massless Dirac operator in the contin- 
uum theory the Atiyah-Singer index theorem Jl5[ 
holds: 



Qtop = index(D) . 
Here 



S2ir z J 



(24) 



(25) 



is the topological charge of the gauge field, and 

index(£>) = n_ — n+ , (26) 

where n + , rt_ is the number of solutions to Du = 
with +/— chiralities, i.e. 75U = ±u. 

The spectrum of D in the continuum consists 

of 

• A = 0; the corresponding modes have defi- 
nite chirality 

• A = ±ia (a real, non-zero); for the corre- 
sponding modes (ua,7sua) = 0. 

Obviously, rt_ — n+ is topological invariant: an 
unpaired eigenvalue at A = cannot be continu- 
ously produced or destroyed. 

In the lattice formulation one has problems 
with both sides of eq. (pi|). On the lhs. there 
are several ways to assign Qtop{U) to a given 
gauge configuration (field theoretical definition, 
geometrical definition). To avoid ambiguities for 
strongly fluctuating gauge fields, cooling is often 
applied. An additional problem is that, in gen- 
eral, dislocations could be present, i.e. configura- 
tions which violate the inequality 



S g (U)>Sf°£-\Q top (U)\. 



(27) 



On the rhs. in an arbitrary gauge field D w has 
no exactly zero eigenvalues - they are distorted 
by lattice artifacts. The spectrum of D w consists 
of 

• real eigenvalues, with (ua,75Ua) ^ but 
not exactly ±1 

• complex eigenvalues, in cc. pairs, A and A* 
with (ma,75"a) = 0. 



For smooth instanton background (p ^ a) D w 
has an isolated real, nearly zero eigenvalue with 
(ua,75"Ua) ~ +1 or —1. However, for a small 
instanton, or a typical MC configuration (say at 
(3 = 6.0) there are a lot of small real eigenval- 
ues. Therefore the definition of topological charge 
based on investigating the spectrum of D w is 
quite ambiguous Pqji7] . 

Let me remind you that the RG approach has a 
theoretically appealing (although difficult in prac- 
tice) way to define the topological charge of a 
given configuration, in spin models [ fl8| and in 
gauge theories Jll| . The solution is the following: 
take the minimizing configuration on the fine lat- 
tice 

C/finc = U(U co &isc) (28) 

and define the FP topological charge on £/coarsc 
through 

Qfop^coarsc) = Qt^° (U Rnc ) . (29) 

It is easy to show that with this definition there 
are no dislocations, i.e. 



S FP (U)>S^-\Q^(U)\ 



(30) 



holds for any configuration. 

As has been shown in Ref. the FP fermion 
operator has nice chiral properties which allow 
one to establish the exact index theorem on the 
lattice. 

Let me discuss the main points. The GW re- 
lation eq. ( |l5| ) (with a = 1) and 75-D75 = 
implies 



D + D^ = DD^ = D^D. 



(31) 



hence D and commute, i.e. D is normal. 
(As a consequence, the eigenvectors correspond- 
ing to different eigenvalues are orthogonal to each 
other.) The spectrum of D (cf. eq. fll8|)) is on the 
circle A = 1 — e lQ . There are three types of eigen- 
values 

• A = 0, with 75 u a = ±ma (denote their num- 
bers by n+ , n,-), 

• A = 2, with 75WA = i"A (with multiplicities 
n' + , n'_), 
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• A = complex, with 75 u\ — u\* and 

(ma,75"a) = 0. 

Note that @ since Tr(7 5 ) = the index of D 
and 2 — D is opposite: 

n , — n_ = Ti- — n + . (32) 
Define the topological charge density by 

q(x) = jtr{<y 6 D{x,x)). (33) 

For the topological charge Q t o P = J2x l( x ) we 
have 

Qtop = iTr(75^) = -^Tr(7 5 (2-^)) 

A 

= n_ — n + = index(D) . 

Here we have used the properties of eigenvectors 
of D listed above. 

Obviously, Qtop is topological invariant since 
the integer n_— n + cannot change smoothly. This 
fact can also be obtained directly from the GW 
relation. 

Equation (|34|) is the index theorem on the lat- 
tice: it states that one can define a (gauge invari- 
ant, real) topological charge density q(x; U) on a 
given gauge configuration which leads to an in- 
teger topological charge. The topological charge 
defined this way equals to the index of D. Both 
sides are defined through D. Is this a tautology? 
No: the important point here is that one has a lo- 
cal topological charge density not just an integer 
number defining the charge. 

It is easy to see that for smooth configurations 

Q(x) = ^ W-tr (F^F p<J ) + 0(a 2 ) (35) 

6ZTT 

since q(x) is a gauge invariant, pseudoscalar, di- 
mension 4 operator and the prefactor is fixed by 
the fact that in the continuum limit a — > one 
should recover the continuum topological charge. 

I would like to stress that any acceptable so- 
lution to the GW relation defines a topological 
charge for which the exact index theorem holds. 
Note, however, that different solutions will give, 
in general, different values for the topological 



charge - they need to agree only for sufficiently 
smooth configurations. From a practical point of 
view the operator D has to provide a topologi- 
cal charge which is robust enough, giving a rea- 
sonable answer also for relatively rough configu- 
rations. To illustrate this point, consider Neu- 
berger's operator in eq. ( pp[) with a small /j,, say 
/i = 0.01. For this choice only very large instan- 
tons are counted as such when defined through 
eq. d3~i| ) - perhaps only those with p/a > 100. 
Smaller instantons "fall through the lattice" since 
the smallest eigenvalue of D w would shift beyond 
fi. Obviously, this would be a bad Dirac opera- 
tor and a bad topological charge - in spite of the 
fact that they satisfy the exact index theorem, 
eq. (0. 

Any acceptable discretization of the Dirac op- 
erator (without requiring to be a solution to the 
GW relation) could be used to determine the cor- 
rect topological charge for sufficiently small lat- 
tice spacing a - the exactly zero eigenvalue of the 
continuum Dirac operator will move only slightly 
away from zero, remaining real, and no other 
nearly zero eigenvalues will appear. The extra re- 
quirement of GW relation forbids the eigenvalue 
to move away from zero for such configurations. 
Hence no extra conditions (beyond acceptability 
and GW relation) are required for the index theo- 
rem to hold. However, as pointed out above, these 
conditions do not guarantee that the topological 
charge defined this way is the 'correct' one - ex- 
cept for sufficiently large instantons. Of course, 
there is no unique way to define the topological 
charge on the lattice, but for not-too-small in- 
stantons (with a radius of a few lattice spacing, 
say) the geometrical definition is natural. For 
practical reasons, one would like to have a GW 
Dirac operator D which reproduces the natural 
topological charge of not too large instantons as 
well. I discuss these points perhaps more than 
necessary, because there is some controversy in 
the literature 

In Ref. Jl(| it has also been shown that for the 
FP Dirac operator one has in addition 



Q gau g cFP = £ ^ ^ RD ^) xx ^ Q ™ (36) 
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i.e. the topological charge defined by D coin- 
cides with the old definition in pure gauge theory. 
In this case, for example, there are no disloca- 
tions, and one expects small lattice artifacts. 

An important point is that D(x,y;U) is dis- 
continuous in the gauge field. The reason is that 
on the lattice one can connect continuously any 
two configurations, but index(D(U)) has a jump 
on the path connecting, say, Qtop = and 1. 

With the FP action the origin of this disconti- 
nuity is that the minimizing configuration U(V) 
becomes ambiguous when the rhs. of eq. (]?]) has 
a degenerate minimum (when the instanton "falls 
through the lattice"). For Neuberger's Dirac 
operator the discontinuity appears when a real 
eigenvalue of D w passes 1. In this case an eigen- 
value of H — 75(1 — D w ) changes sign, and e(H) 
in eq. ( |l9| ) is discontinuous. 

Note that a similar discontinuity appears when 
the fermions are defined in the continuum, in a 
continuous gauge background obtained by inter- 
polating the lattice gauge field J22j. In this case 
the process of interpolation becomes discontinu- 
ous, similarly to the situation with the FP action. 

An important question: does D become non- 
local when this discontinuity occurs? This would 
be a very bad news for the present approach. For- 
tunately, this is not true in general. In the FP ac- 
tion the two solutions U\ and U2 corresponding 
to the degenerate minima are expected to differ 
only locally. Neuberger's D, as shown by Liischer 
23 , 24|] , remains also local when an isolated eigen- 
value of H passes zero, so a discontinuity of D 
does not imply a loss of locality. 

4. The chiral condensate 

Restoring the lattice spacing a, the spectrum 
of D is given by a circle A = - (l — e ia j with 
— 7r < a < ir, crossing the real axis at A = and 
A = 2/a. Correspondingly, the spectrum of D^ 1 
is a straight line parallel to the imaginary axis, 
of the form 1/A = a/2 + ry, —00 < 7 < +00 for 
arbitrary gauge field. 

Denote the unnormalized expectation value of 
0{ip,ip) m a fixed gauge background by 



Then we have 



-Sf(V^) 



0(ip,iP). (37) 



« 4 E(^> F = <i) F El 



(38) 



Obviously, the a/2 part of 1/A is independent of 
the dynamics and has to be subtracted to define 
a proper chiral order parameter. This yields to 



<^x> sub = (^-2iVfiV c 



1 



(39) 



the same as suggested in Ref. ||. 

Adding a mass term m to D would shift the 
whole circle of eigenvalues. Since the real eigen- 
values at A = 2/a also contribute to chiral rela- 
tions it is better to leave them unchanged. This 
is simply achieved by the replacement 



D 



1 



—am j D - in 



(40) 



which not only shifts the circle but also rescales 
its radius appropriately. With this definition of 
m one obtains 



^<l) p = (*(l-iaD],- 



(41) 



This expression is identical to eq. (p9|), i.e. the 
proper order parameter. This equivalent form 
on the rhs. has been suggested first by Chan- 
drasekharan |25|, who obtained it from Liischer's 
exact symmetry. 

5. Exact chiral symmetry on the lattice 

At this moment we are in a peculiar situation: 
all signs of a symmetry are present (Ward iden- 
tities, index theorem, chiral spectrum, . . . ) but 
the action is not chirally invariant. 

Liischer made an important observation that 
in fact the action possesses an exact symmetry. 
Perform an infinitesimal change of variables ip - ► 
tp + ieSip and ip — > V + i^-Stp with a global flavour 
singlet transformation 



! 1 - \al) ; 



Sip 



Sip = ip[l - ^aD ) 75 



(42) 
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From the GW relation eq. (|T^) it follows that the 
action is invariant 



5(ipDip) = . 



(43) 



The corresponding flavour non-singlet transfor- 
mation is given by 



StP = iP{l--aD) l5 T 



(44) 



Now it seems that we have more symmetry 
than possible - the flavour singlet transformation 
should be anomalous. Where is the anomaly hid- 
den? As pointed out by Luscher, the fcrmionic 
integration measure is not invariant under the sin- 
glet transformation in eq. (E3) 



5[dipdil)] = Tv{a lb D)[d^d^] 

= 2JV f (n_ - «+)[##] 



(45) 



That is, the measure breaks the flavour singlet 
chiral symmetry in a topologically non-trivial 
gauge field where index(Z?) = n_ — n + ^ 0. 
The corresponding (global) Ward identity for the 
flavour singlet transformation is 



(<yo) F = 2Jv f i/<o) F , 



(46) 



where v = Q top {U). 

Note that the non-singlet symmetry is not 
anomalous, (SO) F — since eq. ( |45| ) for this case 
contains Tr (a-/ 5 DT) cx tr(T) = 0. 

Using these exact symmetries of the action, the 
Ward identities follow, of course, directly. By tak- 
ing 



(47) 



and a non-singlet transformation one has 
(80) F = i.e. f§ 

lp (l - \aD^j ipj = (m = 0, V finite) .(48) 

This is the order parameter which is expected to 
be broken spontaneously in the thermodynamic 
limit i|: 



lim lim (-0(1- -aD ) tb ) = S ^ . 



(49) 



In Ref. |25[ there is also a relation connected with 
the U(l) problem (massive rj' meson) - this will 
be discussed later. 

Whether the chiral symmetry is indeed broken 
spontaneously in QCD is a delicate dynamical 
question. The chiral invariant Dirac operator pro- 
vides a firm basis to investigate this question by 
avoiding an explicit breaking. 

Note that the Nielsen-Ninomiya no-go theorem 
is avoided here by the non-standard form of the 
lattice chiral symmetry: D is not invariant under 
the naive 75 symmetry, nevertheless it is invariant 
under the generalized transformation. 

6. Chiral decomposition 

In the continuum, chiral symmetry allows one 
to decompose the fermion field into left and right 
chiral components, transforming independently. 
Besides yielding a convenient formalism for vector 
theories like QCD, this decomposition is a neces- 
sary first step towards the chiral gauge theories. 

The matrix 75 anti-commutes with the contin- 
uum Dirac operator and 75 = 1. This allows to 
introduce the chiral projectors 

P± = ^(l±75), (50) 

and the action decomposes as D cont = 
P+D cont P- 



P- Dcont P+ 



On the lattice it is useful to introduce |26 2 
the operator 



7~5 =75(1 - aD) = 75^ , 



(51) 



where D satisfies the GW relation eq. (|l5|). Re- 
member that V is unitary and 75V75 = v^. Then 



-f 5 D = -D% , and % = 1 . 
One can introduce the projectors 

P± = \(l±%) ■ 
They satisfy the relations 
DP+ = P_D , L>P_ = P+D , 
and allow the decomposition 
D = P+DP- + P-DP+ . 



(52) 



(53) 



(54) 



(55) 
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To conclude the definitions let us introduce the 
chiral components of ip: 

i>L = P-i>, ?i = ?P +! (56) 
i> R = P+il>, i> R = i>P-- (57) 

In terms of these components we have 

^Dip = ^ L D^ L + ^ R D^ R , (58) 

analogously to the continuum case. Note that 
eqs. (|56|,|5^) are not symmetric in ip and ip. As 
discussed later, this asymmetric definition seems 
to play an essential role. 

Analogously to the continuum case, the mass 
term should couple the L and R components (but 
not the L,L or R,R components): 

S = Jj l ^j r + Jj R ij L =tf(l- \aD^j $ . (59) 

The corresponding pseudoscalar is given by 

V = $ L $R-$ R 1>L = fysU-\aD\l>. (60) 

The scalar (pseudoscalar) densities S and V will 
often appear in the formulae below. 

Performing a left chiral transformation means 

8ip L = -ip L , 6i/j L =ip L , (61) 
5iP R = , 8$ R = , (62) 

Using this transformation one has 

5V = S, SS = V. (63) 



The combinations in eqs. (p9| , pOD are the nat- 
ural generalizations of the corresponding contin- 
uum expressions. Note that the 1 — \aD fac- 
tor cancels the contribution of the unphysical real 
eigenmode at A = 2/ a. Using these combinations 
of fields, one obtains very simple, compact rela- 
tions, for example: 



£<^>F 



-iVfTr 75 



(64) 



(1 - \aD) 



(l - \am) D 



;i)p = 



■-N f (n+ - n_) (1) F = —N f u (1) F 
m m 



Note that only the A = eigenmodes contribute 
when evaluating the trace. 

It is instructive to consider the flavour singlet 
Ward identity which leads to a relation for the 77' 
propagator: 



77! 

V 



ra V 



1>f.(65) 



If the rf stays massive in the chiral limit ra — ► 
then the lhs. (after averaging over the gauge 
fields) should go to zero. This leads to the inter- 
esting relation [^5) between the topological sus- 
ceptibility and chiral condensate in full QCD 



^ 2 > = -^<s> + o(™ 2 ) 



(66) 



This relation has been obtained previously in the 
continuum p8| , p9| but it is important to have its 
lattice version as well. It shows, for example that 
the topological susceptibility must be a well be- 
haved physical quantity on the lattice. 

For the flavour non-singlet variant of eq. J65| ) 
the term proportional to v 1 is absent since this 
symmetry is not anomalous. Therefore a non- 
vanishing chiral condensate, eq. (^), yields to 
(quasi) Goldstone bosons with mass raj oc raE. 

We conclude this section with a few remarks. 
Note that 75 coincides with e(H) introduced in 
the overlap formalism. Another point worth men- 
tioning is that for the modes with A = one has 
75 = 75 while for the modes at A = 2/ a they 
differ in sign, 75 = —75, as seen from the defini- 
tion eq. (pl[). This observation will be important 
when we consider chiral fermions. 

7. The 6 parameter 

We can now introduce the 9 parameter, as in 
the continuum. Here I shall restrict the discussion 
to some basic points only. 

Define the lattice Lagrangean with a complex 
mass m and the ^-parameter: 

C = C g +ipDip+mip R ip L +m*ip L ip R — i9q{x) ,(67) 

where C s is the gauge action and q(x) is defined 
through D by eq. (33). By changing the variables 



IpL 



(68) 
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one can restore the invariance of the action by 
to — » e M TO and the non-invariance of the measure 
can be compensated by 8 — > 8 — NfCt. Therefore 
the free energy should depend only on the com- 
bination to exp(i#/iVf), as in the continuum. One 
can proceed now e.g. by repeating (now on the 
lattice) the steps of Ref. j28| where e.g. useful 
sum rules have been obtained for the spectrum 
of the Dirac operator. To illustrate this, consider 
the partition function 

Z = Z{8) = J2 z l6vz v = e" F(m ' e) . (69) 

V 

Under some natural assumptions pa] one has for 
the free energy 

F(m, 8) = -FiVfSRe [me i9/Ni ) + 0(m 2 ) . (70) 

Taking then the second derivative with respect to 
8 (for real to) one obtains 

(v 2 ) = VEto^- (71) 

which is equivalent to eq. (p(j). 

As another example, it is easy to show that the 
fermion determinant is given by 

det D{U)= (72) 
(2mrn(V+(l-^) H 2 ) , 

for v > 0, and the product is over the complex 
eigenvalues only. For v < the prefactor is re- 
placed by {2m*)~ u . Eq. (|72| ) differs from its con- 
tinuum counterpart by the presence of the factor 
l-|A| 2 /4. 

8. Cut-off effects with GW Dirac operators 

An important consequence of the exact chiral 
symmetry on the lattice is that there are no 0(a) 
lattice artifacts, to any order of g. The reason is 
the following. The 0(a) lattice artifacts of any 
action could be compensated by adding an extra 
term which is a lattice version of Tpa^F^ip with 
an appropriately chosen coefficient |50[ . The orig- 
inal action in our case has the exact (generalized) 
chiral symmetry of eq. ( [l2"| ) and its artifacts could 



be compensated only by a similarly invariant ex- 
tra term. Any lattice regularization of ipcr^F^ip 
will violate, however, this symmetry, hence the 
0(a) artifacts should be absent for the original 
lattice action as well. Note also that the clover 
term needed for 0(a) improvement has been de- 
termined by requiring that some Ward identity 
holds independently of the gauge field. This is, 
however, automatically satisfied here for m = 
hence there is no need for extra clover term. 

In addition, in the presence of a mass to there 
are also no O(am) artifacts, if the mass is intro- 
duced properly, as in eq. (|67[). Then the symme- 
try to — > —to forbids such term. Indeed, under a 
finite chiral transformation 

ip -* ip' = %i> , "0 -> -iplh (73) 

ij)Dip is invariant while the mass term ^(1 — 
iaZ3)-0 changes sign. Although due to the 
anomaly the partition function is not invariant 
for v ^ 0, this fact is not relevant for our con- 
siderations. Indeed, for m ^ and sufficiently 
large volume one can restrict the system to the 
topologically trivial sector v = j2g] . 

Note, however, that the 0(a 2 ) lattice artifacts 
could still be large at the relevant lattice spacings, 
and could be quite different for different solutions 
of the GW relation. 

Finally, it is interesting to see (although it is 
not necessary after the previous discussion) that 
the absence of 0(a) artifacts on the tree level fol- 
lows directly from the GW relation: expanding D 
in powers of a it implies that 

D = + ac (V + ^a^F^J + 0(a 2 ) , (74) 

(with some constant c) i.e. a term with csw = 1 
is automatically generated. 

9. Properties of Neuberger's Dirac opera- 
tor 

Here we shall discuss some basic properties of 
Neuberger's Dirac operator and its straightfor- 
ward generalizations. 

Consider first the free case. The spectrum E(p) 
of D defined in eq. ( |l4] ) is shown in fig. [|. together 
with the spectrum for the standard Wilson-Dirac 
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Figure 1 . Spectrum of the Wilson-Dirac opera- 
tor (W), and of Neuberger's operator (N). 



Figure 2. Locality of different GW Dirac opera- 
tors for the free fcrmions. 



operator for p = (p, 0,0). The latter is denoted 
by W and given by 



coshi? = 



3 — 2 cosp 
2 — cosp ' 



(75) 



while N denotes the two branches for Neuberger's 
Dirac operator where cosh E(p) is given by 



1 + sin p, and 



3 — 2 cosp 
2(1 — cosp) 



(76) 



(Note that the lower curve extends only up to 
the point where the two branches join each other 
- this peculiar behaviour is due to the non- 
analyticity in the square root.) The spectrum 
for ap > 1 is worse than that for D w . This il- 
lustrates that a solution of GW relation can have 
large 0(a 2 ) artifacts - the fact that D is spread 
over some distance does not necessarily make the 
artifacts small! 

The locality of D for the free case is guaranteed 
by the inequality — y)\\ < C , e~ 7 ' 3: ~ 1 '' where 
\ x ~ v\o = max^ \x — and 7 = 0.693 . . .. 

This asymptotic behaviour is shown on fig. |^ 
by dashed line. The decay of ||D(x)|| for Neu- 
berger's D is shown by circles. Two more solu- 
tions to GW relation are added. Crosses show the 



case with A = 1 — £>hc where -Dhc is a FP ac- 
tion truncated to the hypercube. Since there is a 
non-trivial R in this case the more natural choice, 
A = 1 - VMDhcVZR, (cf. eqs. @,||)) is also 
shown (+ signs). Finally, the true FP action is 
also indicated (triangles). It is obvious that there 
are rather large differences in the range of interac- 
tions for different GW Dirac operators, and this 
is an important factor in making a choice for nu- 
merical simulations (see also Ref. [jL3[). 

A less trivial point is the question of locality in 
a gauge field background. In Ref. |24| it is proven 
that for small gauge fields 



1 — Upi 1 1 < e for all plaquettes 



(77) 



with e < 1/30, the corresponding D(x,y;U) is 
local with exponentially decaying tail, and the 
localization length is uniformly bounded. 

Although eq. (|77]) should hold sufficiently close 
to the continuum limit, it is violated at practi- 
cal values of a. Theoretically one can restrict the 
gauge fields to satisfy eq. (|77|). This is sufficient 
to show that there exists an exactly chiral in- 
variant, acceptable lattice regularization of QCD. 
However, the restriction leads to an astronomical 
correlation length, hence it is only of theoretical 
interest. As mentioned before, in Ref. p3 it is 
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Figure 3. Locality of Ncuberger's Dirac operator 
in gauge field background. 



also shown that an isolated zero of A does not 
spoil locality. 

An important question is the locality in gauge 
fields generated by MC simulations. First results 
on this are presented in Rcfs. [^,^2). In fig. || 
results from Ref. are shown, for SU(3) gauge 
fields generated at (5 = 6.4, 6.2 and 6.0 (triangles, 
squares, circles). The full symbols are for \i = 1 
while the open ones for (roughly) optimized val- 
ues, (j, = 1.2 for (3 = 6.4 and fi = 1.4 for /3 < 6.2. 

It is also important to have a sufficiently ef- 
fective way of evaluating the inverse square root. 
This issue is discussed in refs. p4] , p4] , p2| -p6| . 

I only mention here that it seems to be a fea- 
sible task but still much has to be done in this 
direction. 

A few comments are in order: 

(1) Unlike evaluating D~ x for massless quarks, 
taking the inverse square root is not a critical 
problem - appearance of zero modes of A is acci- 
dental. 

(2) Starting with a good approximation of the 
FP action instead of D w (as in eq. (p2|)) the prob- 
ability to have a zero mode of A is strongly sup- 
pressed, since with the exact FP action A would 
be unitary. Note, however, that for any contin- 



uous parametrization of the FP action, A must 
have a zero mode for specific configurations, when 
an instanton "falls through the lattice". This 
happens on the boundary separating regions of 
gauge field configurations with different topolog- 
ical charges. 

(3) Using a GW Dirac operator there will 
be no "exceptional configurations" in the usual 
sense |37j - these would correspond to eigenval- 
ues A(.D) < 0. (One should distinguish these from 
configurations where X(A) = appearing in this 
approach.) 

(4) One can construct, in principle, configu- 
rations when many "instantons fall through" si- 
multaneously at different sites of a large lattice. 
This would produce a situation when A has many 
nearly zero eigenvalues. Will D be local on such 
configurations? One expects that such configura- 
tions have negligible weight in the partition func- 
tion. Will their effect be also negligible? I think, 
these questions have to be investigated. 

10. Measuring the topological charge 

There are several equivalent ways to measure 
the topological charge defined by D through 
eqs. (@g): 

(1) find the zero eigenmodes of D (or of hermi- 
tian 75 1?), 

(2) count the real eigenvalues of D w in the re- 
gion <= \(D W ) < (i, with chirality signs, 

(3) trace the "level crossings" of H(fi') = 
75 (A*' — D w ) while fi' runs from to /i. 

Note that although we rely here on the exact 
index theorem on the lattice, the value of the 
topological charge determined through D with 
A = /i — D w is the same as the method using 
D w suggested a long time ago |IqJi7|| . The last 
method has also been used in Ref. U in the frame- 
work of the overlap formalism. 

This definition of the topological charge works 
well for small instantons, as expected. But how 
does it work for typical MC configurations? Does 
D (generated using D w ) give a 'better' definition 
than the geometrical one? Hernandez (3£| argues 
that the answer is no. Gattringer and Hipp ]39| 
simulated SU(2) on a 12 4 lattice at /? = 2.4 and 
found that the value of the topological charge, 
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v on individual configurations depends strongly 
on the parameter since there are usually many 
real eigenvalues in the region of interest, and only 
very few configurations have a gap in the spec- 
trum where v does not depend on fj, chosen in 
a reasonable range. They also note that a naive 
inclusion of the clover term makes the situation 
worse. The same conclusion is found by DeGrand, 
A. Hasenfratz and Kovacs p0[ . This last fact is 
not too surprising since the corresponding artifact 
is not an 0(a) effect. 

Note, however, that in the average this defini- 
tion of charge works better: Edwards, Heller and 
Narayanan Q found that the topological suscep- 
tibility stays roughly constant (independent of /i) 
for fi, > 1.0. Here I would like to point out that 
the /^-independence of (v 2 ) does not necessarily 
mean the 'true' value: if v = i^truc + noise then 
(v 2 ) > (^ t 2 ruo ). A good definition of the topologi- 
cal charge (at a given lattice spacing) should give 
relatively unambiguous definition of v on individ- 
ual configurations. In other words, the value of 
the lattice spacing a when a GW Dirac operator 
gives a 'good' definition of the topological charge 
could strongly depend on the choice of D. With 
the presently used lattice spacings the original 
choice of Neuberger (using D w ) seems to produce 
a topological charge which is too ambiguous. 

On the other hand, measuring the topological 
charge with the FP action in the Schwinger model 
po[ gives quite convincing results. For 4d gauge 
theories, however, the situation is less satisfactory 
|l9[ - mostly because of technical problems to 
perform the required multigrid minimization on 
the gauge configurations. 

For QCD a good strategy would be a compro- 
mise: take a Dirac operator Dq with good prop- 
erties (small cut-off effects, nearly satisfying the 
GW relation,. . . ) and build D from it, which in 
turn could be used in the definition of the topolog- 
ical charge. It would have smaller 0(a 2 ) artifacts, 
better locality, fewer 'dangerous' configurations, 
larger gap in the real eigenvalues hence less am- 
biguity of Qtop , but will also need more program- 
ming effort and extra computational overhead. 

Measuring the topological charge density using 
q(x) = ^tr(75D(a;, x; U)) is perhaps too time con- 
suming to be feasible. 



11. On chiral gauge theories 

The question whether it is possible to define 
chiral gauge theories (like the Standard Model) on 
the lattice is a purely theoretical problem, but of 
great importance. Much has been done in this di- 
rection and has been reported earlier |^l]|| . Nev- 
ertheless, there remain several open questions. 

A GW Dirac operator seems to be a good start- 
ing point: it has exact gauge invariance and one 
can define Weyl fermions (cf. eq. ( J5(||57| )). Here I 
shall report briefly about new developments (par- 
tially yet unpublished) by Liischer Q based on 
the properties of GW Dirac operators. The dis- 
cussion will touch only a few interesting points 
which could be easily understood on the basis of 
previous sections. 

Consider fermion fields which are purely left 
handed (see eq. @)) 

P_^ = iP, ^P + =^. (78) 

These are local, gauge invariant conditions. Intro- 
duce a basis for ip(x) and ip{x) through relations 

3 

v k P+ = v k , ip(x) = ^c k v k (x) , (80) 
k 

where Cj and c k are Grassmann variables. In this 
basis the fermion action reads 

S F = i)D^ = c k M k3 c 3 (81) 

kj 

with M kj = ~}Z x v k Dvj. 

The fermion propagator (when there are no 
zero modes of D) is given by 

<V(*My)) F = (1)f (P-D- 1 P + ) . (82) 

It is easy to show that this is chiral in the gener- 
alized sense. Note that 

P-D- 1 P + = P-D- 1 P+ + ^P+ . (83) 

The first term on the rhs. is chiral in the usual 
sense (with 75) while P+ is a contact term. 

Consider a chiral gauge theory with fermion 
multiplet a = 1, 2, . . . , N, with integer charges 
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e Q . For the number of zero modes (n+) Q , (n_) Q 
the index theorem implies 

(n_)„ - (n+) a = e 2 a v (84) 

since FF is replaced by & 2 a FF for smooth gauge 
fields when e a ^ 1. 

A simple but important observation is that 
the number of "0 and ip fields are not necessar- 
ily equal. They are given by Tr(P_) and Tr(P + ), 
respectively. For their difference we have 

Tr (p_) - Tr (P + ) = ^Tr ( 7s £») = „ £ ^ . (85) 

a 

To illustrate this, consider a background field 
with i/ = 1, where n_ = 1, n + = and n'_ = 0, 
n' + = 1. Then for the A = mode 

75" = 75" = (86) 

while for the eigenvector z/ at A = 2 

75"' = -%v' = v (87) 

According to eqs. ( (78^801) , at A = we have only 
ip (but not tp) while at A = 2 both ip and z/> are 
present since 75 and 75 differ by sign at this point. 
Obviously for complex values of A both ip and ip 
are present again. Altogether there is one more 
ip than ip. Of course, this statement depends on 
the topological charge of the gauge field. 

As a consequence the fermionic average (C) F 
of any operator in a given background gauge field 
vanishes unless the fermionic number of type a 
(the number of ip a 's minus the number of ip a 's) 
in O equals to (n_) a - (n + ) a = vJ2 a e i- 

Perform now (for the U(l) case) a global gauge 
transformation g(x) — exp(iu). This modi- 
fies O by a phase factor exp(i(j>) where (p = 
io J2 a e a( n - - n + ) a = w J2 a e a- From the other 
side a global U(l) gauge transformation does not 
change the gauge fields at all, hence invariance 
under this transformation requires 

E e « = > ( 88 ) 

ot 

i.e. the condition of anomaly cancellation of the 
continuum theory. 

These observations, however, constitute only 
the starting point of a difficult mathematical 



problem - to find a satisfactory choice of the 
fermionic integration measure ("chiral determi- 
nant"), i.e. of the basis Vj(x; U). 

Liischer has proven that it is possible to con- 
struct a basis for the U(l) case that satisfies the 
following requirements: 

(1) The fermionic expectation values (0) F are 
smooth functions of U 

(2) The field equations (obtained by some 
change of variables) are local: 

(5C(y)0( Xl )0(x 2 ) . . .) =0 (89) 

if y is far from x±, X2, 

(3) (C) F is gauge invariant if O is. 

Note that the last property may not be neces- 
sary to require. In principle, it is possible that 
after the leading violation of gauge invariance is 
cancelled by the condition e a = 0, the rest 
is small and goes away in the continuum limit, 
without the need of extra counter terms. This 
possibility is based on a simple but surprising re- 
sult of Foerster, Nielsen and Ninomiya |42|| , who 
show that gauge invariance on the lattice is robust 
in the sense that the effect of a sufficiently small 
violating term vanishes in the continuum limit. 
The constructions of the overlap and gauge fixing 
formalisms assume this scenario. It is clear, 
however, that it is much safer - if possible - to 
preserve chiral gauge invariance at finite lattice 
spacing. 

12. Summary 

It has been illustrated that Dirac operators sat- 
isfying the GW relation provide an elegant solu- 
tion to the chirality problem in lattice QCD. In 
particular, there is no additive quark mass renor- 
malization, and Ward identities, soft pion theo- 
rem, etc. are valid, there are no 0(a) artifacts. 

It is pointed out that there is a large freedom 
in choosing a GW Dirac operator, and the choice 
should be optimized for locality, smallness of ar- 
tifacts, etc. 

For practical implementations of these ideas 
new algorithms are needed. That the task is fea- 
sible is illustrated by MC calculations in the do- 
main wall approach H which could be considered 
as a particular approximation to a GW Dirac op- 
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erator. Not restricting itself to the given ideology 
could open ways to more efficient algorithms. 

Finally, as shown by Liischer's construction for 
the U(l) case, it seems to be possible to define 
chiral gauge theories on the lattice based on a 
GW Dirac operator. 
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